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1 Introduction 

\ The CKM paradigm for quark mixing in the standard model accounts for the ob- 



served CP-symmetry violation in the quark sector. The CKM mechanism describes 



D ■ all transitions between quarks in terms of only four parameters: three rotation angles 

and one irreducible phase. Consequently, the flavor sector of the standard model is 
highly predictive. 

^ \ One particularly interesting prediction is that the direct and mixing-induced CP 

I asymmetries, denoted C and S, respectively, are approximately the same in decays 

governed by 6 — )■ qqs (g = u, d, s) transitions and those found in 6 — )■ ccs transitions. 
The level of approximation was estimated as ~ 1% [1]. The dominant contributions 
to 6 — )■ ccs transitions (typically using B^ — J/ipKg) are from tree- level diagrams. 
On the other hand, since flavor changing neutral currents are forbidden at tree level 
in the standard model, the b — )■ qqs transition proceeds almost uniquely via loop 
diagrams (penguins). In these channels, the weak phases in the respective dominant 
contributions are the same as in the standard model, which explains the prediction 



^ , concerning the approximate equality between the CP-violation parameters. 



The CKM angle /3 is related to the mixing-induced CP asymmetry by S = sin(2/3). 
Yet, other amplitudes than penguin may contribute, to some extent, to 6 — )■ qqs 
transitions (e.g., "tree pollution"). The angle extracted from the corresponding S 
parameter is therefore denoted /3eff rather than /3. Penguin diagrams are affected 
by new particles in many extensions of the standard model. As there may be CP- 
violating sources beyond the standard model, these penguin-dominated charmless 
hadronic B decays, are of great interest because they are sensitive to new physics 
effects at large energy scales. 

Various 6 — )■ s dominated charmless hadronic B decays have been studied in order 
to probe this prediction, by comparing the mixing-induced CP asymmetry measured 
for each mode to that measured in 6 — )■ ccs transitions. 
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Here we present recent time-dependent analyses of meson decays to KgKgK^ 
and K^K~K^ from the BABAR collaboration. The former is a CP eigenstate, while 
in the latter the CP content is determined by means of an amplitude analysis. They 
are described in detail in Refs. [2j and [3]. They both benefit from small theoretical 
uncertainties [4], which make the comparison with ccs transitions more meaningful. 

Measurement of 6 — >■ cod transitions such as B^ — > should also yield 

the same value of sin(2/3) as in 6 — )■ ccs, to the extent that the contributions from 
penguin processes may be neglected. The leading and sub-leading order processes 
contributing to — )■ D^*^^ D^*^~ decays are described by a tree diagram and a 
penguin diagram, respectively. The effect of neglecting the penguin amplitude has 
been estimated in models based on factorization and heavy quark symmetry, and the 
corrections are found to be a few percent [5l [6] . Loops involving non-SM particles 
(for example, charged Higgs or SUSY particles) could increase the contribution from 
penguin diagrams and introduce additional phases. A large deviation of the measured 
CP asymmetry S from the value of sin2/3 measured in 6 — )■ ccs transitions or a non- 
zero value of direct CP violation would be strong evidence of new physics. Here 
we present the latest measurement of time-dependent CP asymmetry of partially 
reconstructed B^ — > D*^D*^ decays from the BABAR collaboration [7]. 

The three analyses below use the final BABAR T(4S') dataset that consists of 468 x 
10^ BB decays. 

2 Time-dependent CP asymmetry of KgKgKg 
decays 

The time-dependent CP- violation parameters S and C are extracted by modeling the 
proper-time distribution of B^ — )■ KgKgKg decays with — )■ tt~^7t~, and events 
where one of the Kg mesons decays to tt^tc^. The result is 

S = -0.94 ±0.06, 
C = -0.17 ±0.18 ±0.04, 

where the first quoted uncertainty is statistical and the second is systematic. The 
correlation between S and C is —0.16. 

A two-dimensional statistical likelihood scan in S and C, which is then convolved 
by the systematic uncertainties on the two parameters, is shown in Fig.[Tl We find that 
CP conservation is excluded at 3.8 standard deviations, and thus, for the first time, 
we measure an evidence of CP violation in B'^ — )■ K^KgKg decays. The difference 
between our result and that from B^ — t- ccK^*^ is less than 2 standard deviations. 
The scan also shows that the result is close to the physical boundary, given by the 
constraint S"^ + < 1. 
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Figure 1: Two-dimensional likelihood scan of — 2Aln£ as a function of S and C, 
including systematic uncertainty, in — KgKgKg. The gray scale is given in 
units of a/— 2Aln£. The result of the BABAR analyses of 5° — )■ ccK^*^ decays [8] is 
indicated as a white ellipse and the physical boundary {S^ + < 1) is marked as a 
gray line. The scan appears to be trimmed on the lower left since the PDF becomes 
negative outside the physical region (i.e., the white region does not indicate that the 
scan flattens out at 5cr). 

3 Dalitz-plot analysis of K^K'Kg decays 

Unlike K^K^Kg, K^K~Kg is not a CP eigenstate, and therefore a Dalitz-plot anal- 
ysis is necessary in order to disentangle CP-even from CP-odd states and mea- 
sure /3eff- Because of the importance of understanding the Dalitz-plot structure in 
B^ K+K-K% we study the related modes 5+ ^ K+K-K+ and B+ ^ K^K^K^ 
along with B^ — )■ K~ K^. The mode B^ — )• K^K~K^ is valuable because it has 
the most signal events by far of any B — )■ KKK mode. Far fewer events are observed 
in B^ — K^K^K^ , but its Dalitz-plot has a simplified spin-structure due to the fact 
that the two mesons in the final state are forbidden (by Bose-Einstein statistics) 
to be in an odd angular momentum configuration. A further benefit of a Dalitz-plot 
analysis is that it allows both sin 2/3efr and cos 2/3efT to be determined, through the in- 
terference of odd and even partial waves, which eliminates a trigonometric ambiguity 
between fies and 90° — /3eff- 

We describe the distribution of signal events in the Dalitz plot using an isobar 
approximation, which models the total amplitude as a coherent sum of amplitudes 
from individual decay channels ("isobars"). mesons are reconstructed both in 
K^g TT+TT- and 7r°7r°. 

We vary three sets of /3eff and the CP asymmetry, Acp-, values in the fit: one for the 
0(1020), another for the /o(980), and a third that is shared by all the other charmless 
isobars in order to reduce the number of fit parameters. Note that this last set of 
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isobars contains both even-spin and odd-spin (P-wave) non-resonant terms. Figure [2] 
shows the time- dependent asymmetry for signal- weighted events, both for the 0(1020) 
region (1.01 < rriKK < 1.03GeV/c2) and the 0(lO2O)-excluded region. We find five 
solutions that are essentially degenerate in f3es, and determine f3es{(pKg) = (21±6±2)° 
from a likelihood scan. This is the most accurate result for Pes{4'Kg) at present. 
Excluding the (f){1020)Kg and fo(980)Kg contributions, we measure /3eff = (20.3 ± 
4.3 ± 1.2)° for the remaining — )■ K~^K~Kg decays, and exclude the trigonometric 
reflection 90° — P^s at 4.8cr, including systematic uncertainties. All the results are in 
a good agreement with /3, measured in 6 — ccs decays. 
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Figure 2: The asymmetry {Nbo -N^o)/ (iV^o + N^o ) in 5° ^ K+K 
of At, in the 0(1020) region (left) and 0(lO2O)-excluded region (right). The points 
represent signal-weighted data, and the line is the fit model. 



4 Time-dependent CP asymmetry of partially re- 
constructed B'^ — > D*^D*~ decays 

Since D*^D*~ is the decay of a scalar to two vector mesons, the final state is a 

mixture of CP eigenstates. The CP-odd and CP-even fractions, as well as the corre- 
sponding CP asymmetries have been previously measured from the angular analysis 
of completely reconstructed B^ D*^D*^ decays both by BABAR [9] and Belle [TO] . 
Here, we report a recent measurement from BABAR based on the technique of partial 
reconstruction, which allows one to gain a factor of ~ 5 in the number of selected 
signal events with respect to the most recent BABAR full reconstruction analysis in [9j. 
This result is complementary to the latter measurement, because the statistics used 
are largely independent of each other. 

In the partial reconstruction of a B" — )■ D*~^D*~ candidate, we fully reconstruct 
only one of the two D*^ mesons in the decay chain D* D^n, by identifying 
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candidates in one of four final states: i^vr, Kirir^, KinnT, Kgrnr. Since the kinetic 
energy available in the decay D* — )■ D'^ tt is small, we combine one reconstructed D*^ 
with an oppositely charged low-momentum pion, assumed to originate from the decay 
of the unreconstructed D*^, and evaluate the mass m^ec of the recoiling meson 
and its direction by using the momenta of the two particles. For signal events mrec 
peaks at the nominal mass, while for background events no such peak is visible. 
The direction of the recoiling meson helps to reduce the additional dilution in B 
flavor tagging due to tagging tracks from the unreconstructed D^. 

The analysis proceeds with a series of unbinned maximum-likelihood fits, per- 
formed simultaneously on the on- and off-resonance data samples. The procedure 
can be logically divided in the following three steps: determination of the signal frac- 
tion and several PDF-shape parameters, determination of the tagging dilution due 
to wrong tag assignments, and finally the time- dependent fit to the data, fixing all 
parameter values obtained in the previous steps. 

The final results for C and S, with their correlation coefficient p are: 

C = +0.15 ±0.09 ±0.04 

S = -0.34 ±0.12 ±0.05 ^ u.uu^a. 

The measured values of S and C that we obtain from data only represent a weighted 
average of the CP-even and CP-odd wave function components. If penguin ampli- 
tudes can be neglected then 5*+ = —S-, C+ = — C_ and the value of the CP-even 
components S+ and C+, which we are interested in, can be obtained using the rela- 
tions: 

C = C+ 

S = S+{1-2R^), 

where the factor (1 — 2Rj_) represents the dilution introduced by the CP-odd com- 
ponent R± in the signal. To compute 5+ we use the value measured by BABAR of 
(i?_L = 0.158 ± 0.029) [9j, where the uncertainty is the combined statistical and sys- 
tematic. To evaluate the related systematic uncertainty, we vary this value by ±la. 
We obtain 

C+ = ±0.15 ±0.09 ±0.04 

S+ = -0.49 ±0.18 ±0.07 ±0.04, 

where the uncertainties shown are statistical and systematic; the third uncertainty is 
the contribution from the error on R± described above. This result is compatible with 
previous measurements from BABAR [S] and Belle [ID] using fully reconstructed decays. 
It is approximately 20% more accurate than the previous BABAR measurement. This 
result well agrees with the standard model expectation of negligible contributions to 
the decay amplitude from penguin diagrams and thence with 5+ = — sin 2/3. 



5 



References 

H. -Y. Cheng, C. -K. Chua and A. Soni, Phys. Rev. D 72 (2005) 094003 



[1] 
[2] 

[3] 

[4] 

[5] 
[6] 
[7 
[8 

[9 
[10 



[hep-ph/0506268j . 

J. P. Lees et al. [BABAR Collaboration], Phys. Rev. D 85 (2012) 054023 
|arXiv: 111 1.36361 [hep-ex] ] . 

J. P. Lees et al. [BABAR Collaboration], Phys. Rev. D 85 (2012) 112010 
pEXiv:1201 .5897 [hep-ex]]. 

H. -Y. Cheng, C. -K. Chua and A. Soni, Phys. Rev. D 72 (2005) 014006 
piep-ph/0502235] . 



Z. -z. Xing, Phys. Lett. B 443 (1998) 365 |hep^h/9809496] . 

Z. -z. Xing, Phys. Rev. D 61 (2000) 014010 [ hep-ph/990 7455l. 

J. P. Lees et al. [BABAR Collaboration], l arXiv: 1208. 12821 [hep-ex]. 

B. Aubert et al. [BABAR Collaboration], Phys. Rev. D 79 (2009) 072009 



larXiv:0902.1708 [hep-ex]]. 

B. Aubert et al. [BABAR Collaboration], Phys. Rev. D 79 (2009) 032002 
[ arXiv:0808.18 66 [hep-ex]]. 

K. Vervink et al. [BELLE Collaboration], Phys. Rev. D 80 (2009) 111104 
[arXiv:0901.4057. [hep-ex] ] . 



6 



